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Abstract 



> 

■ Let D C M"^, be a bounded domain with a C°° connected boundary S, L = 
, — V'^ + q{x) be a selfadjoint operator defined in H = Lp'iD) by the Neumann 

boundary condition, ^(x, y, A) be its spectral function, 6'(x, y, A) := ipj{x)ipj{y), 

o : 

■ where Lipj = Xjipj, ^PjnIs = 0, \\^j\\L^(D) = 1) J = 1) 2, . . . . The potential q{x) is a 
^ I real- valued function, q G C°°{D). It is proved that q{x) is uniquely determined by 

Oh' the data 6'(s, s, A) Vs G S, VA G if all the eigenvalues of L are simple. 
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1 Introduction 

Let D C M", be a bounded domain with a connected boundary S", that is, locally 
^ ' the equation x„ = /(x'), x' := (xi . . of S" is given by the function f{x') G C°° . 

■ We assume n = 3. This assumption is used only in the proof of Lemma 3.1 below and 
can be dropped. If > 3 one should refer to the existence of the coordinate system in 
which the metric tensor, used in the proof of lemma 3.1, has the property: (7*" = for 
1 < z < n. The basic ideas of our proof, outlined below, are valid for n > 3. 

Let L = —A + q{x), where A is the Laplacian and q{x) e C°°{D), the potential, is 
a real- valued function. Let fj{x), j = 1,2, . . . , be the normalized real- valued eigen- 
functions of the selfadjoint operator L defined in if = L'^{D) by the Neumann boundary 
condition: 



L(f j = Xjifj in D, \\ipj\\ = l, \\(p\\ := \\ip\\l2(^d), (1-1) 



*key words: Non-overdetermined inverse problem, spectral function, Schrodinger equation 
'I'Math subject classification: 35R30 



LpjN = on 5, 



(1.2) 



where N is the unit normal to S pointing into D' := We choose ipj to be real- 

valued functions. By elliptic regularity, the functions v^j(x) are C°°{D) and C°°{S) on 
the boundary S. 

The spectral function of L is defined by the formula: 

e{x,y,X) := <^j{x)<^j{y), (1.3) 

Aj<A 

the eigenvalues Xj are counted according to their multiplicities. 
The inverse problem (IP) we are interested in is: 
IP: given 9{s, s, A) Ws E S and VA G find q{x). 

We are concerned with the uniqueness of the solution to IP. This inverse problem is 
not overdetermined: the data is a function of n variables s, A ( 5 is an (n — l)-dimensional 
manifold, A G M) and the unknown q{x) is also a function of n variables. 

Most of the non-overdetermined inverse problems are open |Tl[]-|T^, for example: 

a) inverse potential scattering problems with the data A{a' , ao? k) (fixed direction 
of the incident wave, all directions a of the scattered wave and all energies A; > 0), 

b) inverse backscattering problem with the data A{—a, a, fc), 

c) inverse obstacle scattering problem with the data A{a' , ao, fco) (fixed incident di- 
rection, fixed frequency and all directions of the scattered wave), 

and 

d) inverse geophysical problem with the data u{x' , x„ = 0, k) (surface data for all 
frequencies, fixed position of the point source on the surface). 

It seems that no uniqueness theorems for multidimensional inverse scattering-type 
problems with non-overdetermined data have been obtained so far. The IP, studied in 
this paper, is a multidimensional one of the above type and with non-overdetermined 
data. We prove a uniqueness theorem for this problem. 

The inverse problem with the overdetermined data, which is the spectral function 
e{x, y, A), known yx,y e S and VA G M has been considered in [g, [[TT), [|g-[|T^| 



In |Tl||, |T3[, |jTO[ (see also Remark 2.1 below) the following uniqueness theorem is 
obtained: 

Theorem 1.1. The data {\j,{pj{s)}yses determine q{x) uniquely. 

The result we want to prove is: 

Theorem 1.2. The IP has at most one solution if all the eigenvalues of L in (1. !)-(!. 2) 
are simple. 
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It is possible that the conclusion of Theorem 1.2 holds without the assumption con- 
cerning the simplicity of the eigenvalues, but we do not have a proof of this currently. 
Generically one expects the eigenvalues to be simple, if there are no symmetries in the 
problem. 

The strategy of our proof is outhned in the following three steps: 
Step 1. e{s, s, A) =^ {Xj, (/7|(s)}v3e5 

This step is done under the additional assumption that all the eigenvalues Xj are 
simple. 

It is an open question whether Theorem 1.2 holds without this assumption. Any 
assumption that will make Step 1 possible is sufficient for our purpose. 
Step 2. ip]{s) ^ ipj{s) 

Step 3. Apply Theorem 1.1 to the data {Xj, i^j{s)}ys^s 

Step 1 does not require much work. 
Step 2 requires the basic work. 

Let us outline the ideas needed for the completion of Step 2. 

One can prove that ^Pj{s) cannot have zeros of infinite order, that is, if y G S* and 
< Cm\s — y\™' Vm = 1,2, ... , where Cm are positive constants and s & S are 
arbitrary points on S, then (fj{x) = 0. In the Appendix it is proved that if a real- valued 
function / G C°°(5') does not have zeros of infinite order on S, then, up to its sign, the 
function / is uniquely defined on S by its square p. A procedure for finding, up to its 
sign, ^Pj{s) on S from the knowledge of V'j(s)^ on S is described below. 

It is known that ^j{s) cannot vanish on an open (in S) subset of S (uniqueness of the 
solution to the Cauchy problem for elliptic equations). Thus, '^j{s) ^ for almost all 
points of S. Take an arbitrary point sq at which </^j(so) 7^ 0. By the continuity of V'j(s) 
on S, there is a maximal domain T <Z S containing the point sq in which (p'j{s) ^ 0. 

Define (fij{s) :— y^^^(s) > in J^. Let C (Z S he the boundary of J^. We prove that 

(fj{s) is uniquely determined in J-'' := S\J-. The problem is to determine the sign of 
(Pj{s) for s G JF' in a neighborhood of C, where fj{s) = when s G C 

The basic idea for determining this sign is to extend (fj{s) = y''^|(s) from T across 
C with maximal smoothness. Such an extension is unique and determines the sign of 
(Pj{s) outside of JF. For example, if one defines y/x^ = x > ioi x > 0, then the unique 
maximally smooth extension of Vx^ in the region x < is Vx^ = x. 

To determine the sign of (pj{s) outside of J^, calculate the minimal integer m for 

Vji^) — y^<y?j(s), s e !F G S, such that 7(5) := lims^s-00: 7^ 0. Here s ^ s e JC 
along a curve £ originating at a point in transversal to JC, and passing through a point 
s G £ into J^'. 

The integer m < cxd does exist if the zeros of fj{s) are of finite order. That these 
zeros are indeed of finite order is proved in Lemma 3.1 in section 3. 

Let us describe the way to continue (pj{s) along the curve i across C into J^'. 
Define sgn (pj{s) — (—1)"*, where s G C is any point in a sufficiently small 
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neighborhood of s. This algorithm determines uniquely ifj{s) for all s & S, given the 
data v^^(s). Since the eigenvalues are assumed simple, the above algorithm produces the 
trace (Pj{s) of the eigenfunction (pj{x) on S. 

Since q{x) and 5* are C°°-smooth up to the boundary, so are the eigenfunctions (Pj{x), 
and the above algorithm produces a C°°-smooth function on 5*. Therefore one can use 
the Malgrange preparation theorem (0, p. 43) to study the set of singular points of C. 

The above argument deals with the continuation of an eigenfunction through its "zero 
line" C on 5*. There are at most finitely many points on a compact surface S at which 
several "zero lines" intersect and only finitely many "zero lines" can intersect at one 
point. Otherwise there would be a point on S which is a zero of infinite order of (fj{s), 
and this is impossible by Lemma 3.1 of section 3 below. 

If one continues v^j(s) across C by the above rule, the function fj{s) will be uniquely 
determined on all of S by the choice of its sign at the initial point sq. 

In the Appendix, at the end of the paper, we include an alternative proof of a state- 
ment we used in a discussion of Step 2. It is proved in the Appendix that a smooth real- 
valued function, which is defined on a smooth connected manifold M without boundary, 
and has no zeros of infinite order on M, is uniquely, up to a sign, determined on M by its 
square. This proof was communicated to me on behalf of Professor Yu. M. Berezansky. 

This completes the description of Step 2. 

The C°° smoothness of the data is assumed for technical reasons: it guarantees the 
existence of an integer m in the above construction. It would be interesting to weaken 
this assumption and to find out if one can prove Lemma 3.1 assuming S is C^'^— smooth 
and q E 

In section 2 we sketch for convenience of the reader a proof of Theorem 1.1. In section 
3 a detailed discussion of Step 2 is given. 



2 Proof of Theorem 1.1 

The argument of this section is based on the ideas from (||13||) 

Without loss of generality assume that Xj ^ Vj and take in what follows n = 3. 
The argument for n > 3 is essentially the same. Let 

Lu := —Au + q{x)u = in D, ujy = f on S. (2.1) 

Then 

^(^) = ^ y (2-2) 

j=i 1 

Assume that ^1(0;) and q2{x) generate the same data 

{Aj-,^j(s)}vj. (2.3) 
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Take an arbitrary / G iJ2 (5'), denote p{x) := q2{x) — qi{x), w{x) := Ui{x) — U2{x), 
where Um{x) is the solution to (2.1) with q = qm{x), and subtract from equation (1.1) 
with q = qi equation (1.1) with q = q2. The result is 

Liw := —Aw + qiw = pu2, on S, (2.4) 

where Lm '■= — + qm{x), m = 1,2. Clearly wn = on S*. We will prove later that 

w = on 5. (2.5) 

Assuming (2.5), multiply (2.4) by an arbitrary vi G N{Li) := {vi : Liv = in f i G 
H'^{D)}, where H'^{D) is the Sobolev space, integrate by parts, use the boundary condi- 
tions (2.4) and (2.5) and get the orthogonality relation: 

= I p{x)u2{x)vi{x)dx WvieN{Li), Wu2eN{L2). (2.6) 
Jd 

By property C ( see |11[), (2.6) implies p{x) = 0, i.e. gi(x) = q2{x). 
Recall that property C means that the set of products {f i(a;)-U2(a^)}vt)ieAf(Li)VM2GAf(-L2) 
is complete in U'{D),p > 1. A detailed discussion and proof of this property is given in 

ini- 

To complete the proof of Theorem 1.1, note that (2.5) follows from (2.2) and from 
the basic assumption that qi and q2 generate the same data (2.3). From (2.2) it fol- 
lows formally that Ui{s) = U2{s). Indeed, if ipf\s) = ^f\s) and A^"^'' = A^-^"* Vj, 

where y^j"*^ and A^'"'' are the eigenf unctions and eigenvalues of the operator Lm,Tn = 
1, 2, with q = qm, then formula (2.2) with x = s & S shows that ui = U2 on S. 

However, the series (2.2) converges in L'^{D) but not necessarily on S. Therefore a 
justification of the above argument is needed. Such a justification can be given folllowing 
an idea in [10]. Theorem 1.1 is proved. □ 

3 Finding (Pj{s) from (p'j(s). 

In the introduction a method for finding (pj{s) from the knowledge of v^j(s) has been 
discussed. 

Here a justification of this method is presented. This justification consists mainly of 
the proof of the following lemma. 

Lemma 3.1. The function ipj{s) does not have zeros of infinite order. 

Proof. It is known (see |§, p. 14, where a stronger result is obtained) that a solution to 
the second order elliptic inequality: 

\Mou\ < c{\u\ + \gradu\) in D, (3.1) 
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where c = const > 0, and Mq is a strictly elliptic homogeneous second order differential 
expression (summation is understood over the repeated indices): 



MoU := -gmjUmj, Umj ■= , (3.2) 

cannot have a zero of infinite order at a point y E D provided that: 

gmj{x) are Lipschitz-continuous and strictly elliptic, (3.3) 

and 

gmj{0) are real numbers. (3.4) 
By zero of infinite order of a solution u of (3.1) a point y ^ D is meant such that 

/ \u{x)fdx < c^e"" Ve > 0,Vm = 1,2,...., (3.5) 

J\y-x\<e 

where Cm are positive constants independent of e, and e > is sufficiently small so that 
the ball B{y, e) := {x : \x — y\ < e} G D. 

We use this result to prove that the same is true if S* G C^'^ and y E S. Let x„ = 
f{x'),x' = (xi, . . . , Xn-i), be the local equation of S" in a neighborhood of the point y E S, 
which we choose as the origin, and be directed along the normal to S pointed 

into D' := W\D. 

Let us introduce the new orthogonal coordinates 

ij = ij{x^}, l<j,i<n, 

so that ^„ = is the equation of S* in a neighborhood of the origin, and assume n = 3. 
For example, one can use the coordinate system in which the z— axis is directed along 
the outer normal to S, and the x, y— coordinates are isothermal coordinates on S, which 
are known to exist for two-dimensional C°°— surfaces in ( and even for C^— surfaces 
in R^, [Q, p. 246). For an arbitrary C°°— surface in one can prove that locally one 
can introduce (non-uniquely) the coordinates in which the metric tensor is diagonal in a 
neighborhood of S. To do this, one takes two arbitrary linearly independent vector fields 
tangent to S and orthogonalize them with respect to the Euclidean metric in R^ using 
the Gram-Schmidt procedure, which is always possible. The third axis of the coordinate 
system, which we are constructing, is directed along the normal to 5* at each point of 
the patch on S. Let A and B be the resulting orthogonal vector fields tangent to S. 
Then one can find (non-uniquely) a function x, defined on the local chart, such that the 
Lie bracket of the vector fields wA and B vanishes, [xA, B] = 0. Then the flows of the 
vector fields and B commute and provide the desired orthogonal coordinate system 
in which the metric tensor is diagonal (0). The condition [xA, B] = can be written as 
the following linear partial differential equation: x^x — hxx — = 0, where h := —j^^, 
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X and y are the parameters, and z = f{x,y) is the local equation of the surface S on 
the chart. The above equation for x has many solutions. One can find a solution x{x, y) 
by the standard method of characteristics. A unique solution is specified by prescribing 
some Cauchy data, which geometrically means that a noncharacteristic curve through 
which the surface x = xi^^v) passes, should be specified (@). 

If n > 3 the situation is less simple if one wants to use the same idea in the argument: 
there are many {{n — l){n ~ 2)/2) Lie brackets to vanish in the case of — 1 vector 
fields tangent to S in R^, and it is not clear for what S these conditions can be satisfied. 
However, for our argument it is sufficient to have the coordinate system in which the 
metric tensor has zero elements g^"' for 1 < j < n — 1, and g"'"' does not vanish. In 
this case the even continuation (3.8)-(3.10), that is used below, still allows one to claim 
that the function (3.8) solves the same equation in the region ^„ > as it solves in the 
original region ^„ < 0. If the elements g^" for 1 < j < n — 1 do not vanish, then the 
equation in the region ^„ > will have some of the coefficients in front of the second 
mixed derivatives g^^Wjn with the minus sign, while these coefficients in the region ^„ < 
enter with the plus sign. So, in this case the principal part Mq of the operator M, which 
is used in (3.11), will be different in the regions ^„ > and ^„ < 0. Recall that the 
Laplace operator in the new coordinates has the form of the Laplace-Beltrami operator: 
Aw = g~^^'^{g^^^g^-'wj)i, where g := det{gij), summation is understood over the repeated 
indices, and Wj := 

If y & S is a zero of <fj{s) of infinite order in the sense |v5j(s)| < Cm\s — y\"^ Vm = 
l,2,....,s G S, then equations (1.1) and (1.2) imply that |v5j(a;)| < Cm\x — y\"^ Vm = 
1,2,...., a; G -D n B{y,e), so that condition similar to (3.5) holds for the integrals over 
D n B{y,e). Indeed, the derivatives of (pj{x) in the tangential to S directions at the 
point y & S vanish by the assumption. From (1.1) and (1.2) it follows that the normal 
derivatives of <fj{x) of the first and second order vanish at y. Differentiating equation 
(1.1) along the normal one concludes that all the normal derivatives of fj{x) vanish at 
the point y. Thus we may assume that ?/ G S" is a zero of (pj{x) of infinite order, that is, 
inequahties (3.5) hold for u = (pj{x) for the integral over D fl B{y, e). 

In the ^- coordinates one writes the equation for (fj-. 

—Aipj + q{x)ipj — Xjipj{x) = in D, ipj^ = on S*, (3.6) 

as follows: 

M^j = for ^„ < in D, (^^^^ = at ^„ = 0. (3.7) 

We drop the subscript j of (pj and of \j in what follows. In (3.7) the operator M is 
defined as: 
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over the repeated indices summation is understood, gij is the metric tensor of the new 
coordinate system, g := det{gij), and the coefficients in front of the second-order deriva- 
tives in M are extended to the region ^„ > as even functions of C,n, so that the extended 
coefficients are Lipschitz in a ball centered a.t y & S with radius e > 0. 
Let us define w in a. neighborhood of the origin, |^| < e, by setting 

--(t/^l '!?"-n' (3-8) 
1 , if a > 0, ^ ^ 



and 



The functions gjjiO^ defined by (3.10), are Lipschitz in the ball 5(0, e), if gjj{^) is 
Lipschitz in B{0, e) n i?!!, where R"^ := : ^ < 0}. 
Furthermore, 

\Mow\ <c{\w\ + \gradw\) in 5(0, e), (3.11) 

where c = const > 0, and 

\w\^d^<Cme'^ Vm = l,2, , (3.12) 

JB{0,e) 

if one assumes 

\ip\'^d^ <Cme"' Vm=l,2, , (3.13) 



B(0,e)nlR!i 

for all sufficiently small e > 0. 

Note that the change of variables x — > is a smooth diffeomorphism in a neighborhood 
of the origin, which maps the region Dr\B{0,e) onto a neighborhood of the origin in in 
^— space, and one can always choose a half-ball -8(0, e) belonging to this neighborhood, 
so that (3.13) follows from (3.5). 

Therefore w = in i?(0, e) and consequently (f = in B{0,e) flM" . This implies that 
if = in D hj the unique continuation theorem (see ^ and |T^). This is a contradiction 
since \\(p\ \ = 1. 

Lemma 3.1 is proved. □ 

This lemma provides a justification of the argument given for Step 2 in the introduc- 
tion. 
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Remark 3.1. In this remark we comment on the numerical recovery of the potential from 
the data 9{s,s,\). We have explained how to get the data {\j,(pj{s)}\/j from 9{s,s,\). 
Therefore the spectral function 

e{s,t,X) = ^A^Wt), ^s,teS (3.14) 

Aj<A 

and the resolvent kernel 

G{s,t,X)= rMi^ (3.15) 

are known on S. IfG{s,t, A) is known, then the Neumann-to-Dirichlet (N-D) map un- : = 
f —>■ h := u\s is known. This map A : f h at every fixed X ^ Xj Vj, associates to 
every f G if 2 (5*) a function h G if 2 (5) by the formula 



h{s) = / G{s,t,iJ,)f{t)dt. (3.16) 
Js 

Here 

Lu — Xu = in D, = f on S. (3-17) 
The unique solution to (3.17) is given by (3.16). 

We now outline an alternative proof of Theorem 1.2: steps 1 and 2 are the same as in 
Theorem 1.2; step 3 consists of the construction of the N-D map by formula (3.16) and 
reference to Theorem 3.1 below. 

In [0] and the result similar to Theorem 3.1 (see below) is proved for the D-N 
(Dirichlet-to- Neumann) map. Our proof follows the idea of the proof in [14|. As above, 
Lm, m = 1, 2, stands for the Neumann operator — A + g^ in D. Let us assume that X E R 
is a regular point for Li and L2. 

Theorem 3.1. If Li — X and L2 — X generate the same N-D map then qi = q2. 
Proof. Take an arbitrary / G H^{S) and consider the problems 

{Lm - X)Um = 0m D, U.mN = / OU 5*. 

Subtract from the equation with m = 1 the equation with m = 2 and get 

{Ljn — X)w = pu2 in D, wjy = w = on S, (3.18) 

where w := Ui — U2,p '■= q2 — qi- The condition w = on S follows from the basic 
assumption, namely the N-D map is the same for Li and L2. Let ijji G N{Li) := {ip : 
Liip = in D,?/' G H^{D)} be arbitrary. Multiply (3.18) by ipi, integrate by parts using 
boundary conditions (3.18), and get the orthogonality relation: 

pu2^idx = VV-i G N{Li) Wu2 G N{L2). (3.19) 

D 

The last inclusion in (3.19) follows since / G H^{S) is arbitrary. From (3.19) and 

property C (see |TT[|) it follows that p{x) = 0. Theorem 3.1 is proved. □ 
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In |TT] |T^-|T^] one finds a discussion of the inversion algorithm. 
Appendix 

Lemma. Let M be a connected C°° manifold and suppose that a real-valued function 
f G C°°{M) does not have zeros of infinite order, i.e., for an arbitrary point xq G M, 
f^^K^o) 7^ for some a. Then, up to sign, the function f is uniquely defined on M by 
its square f^. 

Proof. Let g G C°°{M) be a real- valued function satisfying the condition g'^ = 
Then the lemma states that either g = f or g = —f. Consider the function u = f — g. If 
M = on M, then g = f and the lemma is proved. Consider X = {x : u{x) 7^ 0, x G M}. 
This is an open set and u{x) = 2f{x) on this set, since, for a fixed x, the condition 
g'^{x) = f'^{x) implies that f{x) = g{x) or g{x) = — /(x). The equality u{x) = 2/(x) 
will also hold on the closure X. If X = M, then u = f — g = 2/, i.e., g = — / on M, 
and the lemma is proved. Let Y = M \ X he a nonempty open set. Since M = X U Y 
and M is connected, there exists a point Xq E X HY . We have u{x) = 2f{x) on X, and 
the function / does not have a zero of infinite order, hence there exists an a such that 
u^"\xo) = 2f^"\xo) 7^ 0. On the other hand, Xq E Y and u{x) = on the open set Y. 
Thus u'^^\y) = for all y E Y. This and the continuity of imply, as y — xq, that 
= 0. This contradiction proves the lemma. □ 
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